Local quantum criticality in confined fermions on optical lattices 
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Using quantum Monte Carlo simulations, we show that the one- dimensional fermionic Hubbard 
model in a harmonic potential displays quantum critical behavior at the boundaries of a Mott- 
insulating region. A local compressibility defined to characterize the Mott-insulating phase has a 
non-trivial critical exponent. Both the local compressibility and the variance of the local density 
show universality with respect to the confining potential. We determine a generic phase diagram, 
that allows the prediction of the phases to be observed in experiments with ultracold fermionic 
atoms trapped on optical lattices. 
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The Mott metal-insulator transition (MMIT), a 
paradigm of strong correlations, was recently realized in 
ultracold atoms confined on an optical lattice 0. Due 
to the fact that the atoms interact only via a contact 
potential, this system constitutes the most direct exper- 
imental realization of the Hubbard model which is the 
prototype generally used to study the MMIT. Whereas 
optical lattices contain bosonic atoms, recent progress in 
cooling techniques allow fermionic systems to go well be- 
low the degeneracy temperature 1^ y| , such that even su- 
perfluidity appears within reach It is therefore, to be 
expected that soon a fermionic MMIT will be realized on 
an optical lattice, offering the possibility to confront in a 
controlled way our knowledge of the MMIT in solid-state 
systems, without extrinsic effects always present there. 
This possibility is especially important since the MMIT 
is not only a long standing problem in condensed matter 
physics, but has also received renewed attention in recent 
years due to its different manifestations in a number of 
transition metal oxides, the most prominent being high 
temperature superconductors 0. 

Motivated by the possibility of such crossfertilization, 
we performed QMC simulations for the ground-state of a 
one-dimensional Hubbard model with a harmonic poten- 
tial, as in experiments with ultracold atoms, confining 
spin 1/2 fermions. The one-dimensional case was cho- 
sen since in one dimension, the quantum critical prop- 
erties for the unconfined system are well characterized 
by the Bethe Ansatz solution where, in particular, the 
global compressibility k ~ dn/dfi diverges as d~^, where 
S = 1 — n, and n is the expectation value of the density 
However, as shown theoretically and numerically 
in the presence of a confining potential the Mott- 
insulating phase is restricted to a domain that coexists 
with a compressible phase. This is in contrast to the 
global character typical of solid state systems. We show 



in this Letter, that a properly defined local compressibil- 
ity displays critical behavior on approaching the edges of 
the Mott-insulating phase, revealing a new critical expo- 
nent. Furthermore, it is shown that both the variance of 
the local density, =< nf > — < rii >^, and the local 
compressibility as functions of the local density n^, are 
independent of the confining potential for — ^ 1. The 
exponents are also universal with respect to the strength 
of the interaction. 

The Hamiltonian studied is as follows: 

H = -t^^ (clffCj+ia + h.c}j + U ^riiiriii 

where c^^ and Cj^. are creation and annihilation operators, 
respectively, for a fermion on site i with spin <t =t, i- 
The local density per spin is Uia = cl^^Cj^.. The contact 
interaction is repulsive (C/ > 0) and the last term mod- 
els the potential of the magneto-optic trap. The QMC 
simulations were performed using a projector algorithm 
1^ Hol ITll IT ^ . which applies exp{—9H) to a trial wave- 
function (in our case the solution for ?7 = 0). A projector 
parameter 6 ~ 20 /t suffices to reach well converged val- 
ues of the observables discussed here. A time slice of 
At = 0.05/t was used in general. 

Figure n shows density profiles along a harmonic trap 
{a = 2) for different fillings such that the system goes 
from an entirely metallic phase to a phase with insulat- 
ing regions due to full occupancy of the sites, coexisting 
with metallic regions. At this point, such an identifica- 
tion is only based on the occupation number, using our 
knowledge from unconfined periodic systems. Although 
it is natural to identify the regions with n = 1 as insulat- 
ing phases, the global compressibility used in unconfined 
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systems is not a useful order parameter to characterize 
the phases due to the coexistence of local compressible 
regions with incompressible ones. A first quantity that 
can be used instead is the variance of the site density 
Ai, since on entering the Mott-insulating region a sup- 
pression of double occupancy should occur, leading to a 
decrease of the variance. 




FIG. 1: Density profiles along the trap for different fillings. 
Flat terraces are the Mott insulating regions. 

Figure|21(a) shows three characteristic density profiles. 
In all of them, U/t = 6 but depending on the filling (Nf) 
and strength of the potential V2, we obtain (i) an ap- 
proximately parabolic density profile, indicating that the 
whole system is in a metallic phase, (ii) Increasing the 
number of particles, a Mott-plateau develops in the cen- 
ter of the system, and finally, (iii) with still a higher 
filling, a new metallic phase develops in the center of the 
plateau. The potential V2 was varied in order to obtain 
well developed phases. The corresponding profiles for 
the variance A are shown in Fig[2Ib). In general, a sup- 
pression of A is present in the regions where the density 
profile shows a plateau. However, although in the region 
with Ui = \ the variance is lower than in the regions 
surrounding it, it does not vanish. Therefore, while the 
decrease of A is a signature for the Mott-insulating phase, 
still a clearer distinction is needed. For this purpose we 
introduce a local compressibility defined as follows: 

\j\<KU) 

where Xi.j is the density-density correlation function. We 
take the length scale 1{U) ~ b^{U), where £,{U) is the 
correlation length given by Xi.j in the unconfined system 
at half filling for the given value of U. The value of b 
is such that becomes insensitive to the value chosen. 
In general we have & ~ 5 — fO, with (,{U) ~ a (a is the 
lattice constant) for the values of U used here. The lo- 
cal compressibility thus gives the response to a constant 
shift of the potential over a finite range but over distances 
larger than £,(U) in the periodic, unconfined system. If 
a region is in a Mott-insulating phase, and hence incom- 
pressible, no density response over distances larger than 
^ is expected, leading to k\ = 0. Figure |2Ic) shows the 



profile along the trap of k', where the compressibility 
becomes zero in the outlying regions, where no particles 
are present and also where a Mott plateau is present. 
Therefore, the local compressibility defined here serves 
as a genuine local order parameter to describe the insu- 
lating regions that coexist, in general, with a surrounding 
metallic zone or even with metallic intrusions, beyond the 
intuitive pictures on the basis of the density profiles. 




FIG. 2: (a) Density profiles for (A) Nf = 30, f/ = 6t and 
V2 = 15t, (O) Nf = 70, U = 6t and V2 = 6.25t, and (v) 
Nf = 74, f/ = 6t and V2 = 7t. (b) Variance of the local 
density (c) Local compressibility as defined in eq. (|5J. 

Now that phases can be characterized quantitatively, 
we concentrate on the regions where the system goes from 
one phase to another. Criticality can arise, despite the 
microscopic spatial size, due to the extension in imagi- 
nary time that reaches a thermodynamic limit at T = 0, 
very much like the case of the single impurity Kondo 
problem |l3j , where long-range interactions in imaginary 
time appear for the local degree of freedom as a result 
of the interaction with the rest of the system. Recent 
experiments leading to a MMIT consider a system 
with linear dimension ~ 65a, i.e. still in this microscopic 
range. An intriguing future question, for both theory and 
experiment, will be the role of spatial dimension in the 
critical behavior of systems in the thermodynamic limit. 
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FigurcOshows the local compressibility vs. S for (5 — > 
in a double logarithmic plot. A power law ~ is 
obtained, with tu < 1, such that a divergence results 
in its derivative with respect to n, showing that critical 
fluctuations are present in this region. Since the QMC 
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FIG. 3: The local compressibility k' vs. S — 1 — n at S for 
(A) Nf =70, U = 8t and V2 = 6.25t; (v) Nf = 70, U = 6t 
and V2 = G.25t; (Q) Nf =72, U = 6t and a quartic potential 
with V4 = 6.25t; (0) unconfined periodic system with U = 6t. 
Inset: Dependence of the critical exponent vj on Ar^. 

simulation is affected by systematic errors due to dis- 
cretization in imaginary time, it is important to consider 
the limit At in determining the critical exponent. 
The inset in Fig. |21 shows such an extrapolation lead- 
ing to ti7 ~ 0.68 — 0.78. At this point we should re- 
mark that the presence of the harmonic potential allows 
the determination of the density dependence of various 
quantities with unprecedented detail on feasible system 
sizes as opposed to unconfined periodic systems, where 
systems with 10'^ — lO'' sites would be necessary to allow 
for similar variations in density. In addition to the power 
law behavior. Fig. Oshows that for 5^0, the local com- 
pressibility of systems with a harmonic potential but dif- 
ferent strengths of the interaction or even with a quartic 
confining potential, collapse on the same curve. Hence, 
universal behavior as expected for critical phenomena is 
observed also in this case. This fact is particularly im- 
portant with regard to experiments, since it implies that 
the observation of criticality should be possible for real- 
istic confining potentials, and not only restricted to per- 
fect harmonic ones, as usually used in theoretical calcu- 
lations. However, Fig. |3| shows also that the unconfined 
case departs from all the others. Up to the largest sys- 
tems we simulated (400 sites), we observe an increasing 
slope rather than the power law of the confined systems. 

Having shown that the local compressibility displays 
universality on approaching a Mott-insulating region, we 
consider the variance A as a function of the density n 
for various values of U and different confining potentials. 
Figure^lshows A vs. n for a variety of systems, where not 
only the number of particles and the size of the system are 
changed, but also different forms of the confining poten- 
tial were used. Here we considered a harmonic potential, 
a quartic one, and a superposition of a harmonic, a cubic 
and a quartic one, such that even reflection symmetry 



across the center of the system is broken. It appears at 
first glance that the data can only be distinguished by the 
strength of the interaction U , showing that the variance 
is rather insensitive to the form of the potential. The 
different insets, however, show that a close examination 
leads to the conclusion that only near n = 1 and only 
in the situations where at n = 1 a Mott-insulator ex- 
ists, universality sets in. The inset for n around 0.6 and 




FIG. 4: a) Variance A vs. n for (O) harmonic potential 
V2 = G.25t with N = 100; (A) quartic potential V4 = 15.82* 
with A*' = 150; (v) harmonic potential V2 = lOt -\- cubic 
V3 = 2.5* + quartic V4 = 7.5* with = 50; and (full line) 
unconfined periodic potential with N = 102 sites. The curves 
correspond from top to bottom to U/t = 0, 2, 4, 8. For a 
discussion of the insets, see text. 

U = 8*:, shows that the unconfined system has different 
variance from the others albeit very close on a raw scale. 
This difference is well beyond the error bars. Also the in- 
set around n = 1 and for U = At, shows that systems that 
do not form a Mott insulating phase in spite of reaching a 
density n = 1, have a different variance from those having 
a Mott-insulator. The features above show that even a 
very local quantity like the variance cannot be accurately 
described using a local density approximation (Thomas- 
Fermi) [ll, and can lead to even qualitatively wrong 
results, as for U = 4:t and n = I, where such an approxi- 
mation would predict a Mott-insulator instead of a metal 
as in our simulations. Only the case where all systems 
have a Mott-insulating phase at n = 1 (J7 = 8*), shows 
universal behavior independent of the potential, a univer- 
sality that encompasses also the unconfined systems. For 
the unconfined system, the behavior of the variance can 
be examined with Bethc-^nsate (l5i | in the limit 5—^0. 
In this limit and to leading order in 6, the ground state 
energy is given by 0| Eq{5)/N -Eq{5 = 0)/7V cx 5, such 
that the double occupancy, which can be obtained as the 
derivative of the ground-state energy with respect to U, 
will also converge as 5 towards its value at half-filling. 
Such behavior is also obtained in our case as shown by 
the inset sA n = \ {U — 8t) in Fig. 0] Detailed data for 
the variance close to n = 1 will be presented elsewhere. 
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Finally we consider the phase diagram of the system. 
As shown in Fig. |21 unlike the exponents, phase bound- 
aries seem to be rather sensitive to the choice of potential, 
number of particles and strength of the interaction. As in 
the unconfined case, we would expect to be able to relate 
systems with different number of particles and/or sizes 
by their density. Given the harmonic potential, a char- 
acteristic length (in units of the lattice constant) is given 

by N (4V2 /t) , such that a characteristic density can 
be defined. Figure [31 shows that the characteristic den- 
sity p = p\/4:V2/t {p = Nf/N) is a meaningful quantity 
to characterize the phase diagram. There, the phase dia- 
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FIG. 5: Phase diagram for a system with A'^ = 100 (v) and 
A*' = 150 (i i) sites. The phases are explained in the text. 

grams for two systems with different sizes {N — 100 and 
= 150) and different strength of the harmonic poten- 
tial (V2 = 15t and V2 = 11.25t respectively) are depicted 
showing that such a scaling allows to compare systems 
with different sizes, different number of particles, and 
different strength of the potential. This makes it possi- 
ble to relate the results of numerical simulations to much 
larger experimental systems. The different phases ob- 
tained are: A pure metal without insulating regions (A), 
a Mott-insulator at the center of the trap (B), a metallic 
intrusion at the center of a Mott-insulator (C), a "band 
insulator" (i.e. with n = 2) at the center of the trap sur- 
rounded by a metal (D), and finally a "band insulator" 
surrounded by a metal, surrounded by a Mott-insulator 
with the outermost region being again a metal (E). Two 
features are remarkable here. The first one is that on 
varying the filling of the trap, a reentrant behavior is ob- 



served for the phase A. The density profile shows a shoul- 
der as can be seen in Fig. ^ before reaching the plateau 
with n ~ 2 but, as shown by the inset of Fig. ^for U ~ 4 
around n = 1, it is possible to go through a region with 
n = 1 without reaching the value of the variance that 
corresponds to a Mott-insulator. The second intriguing 
feature is that the boundary between the regions A and 
B remains at the same value of the characteristic density 
for all values of U that could be simulated. 

In summary, on the basis of QMC simulations of the 
Hubbard model with a harmonic potential, we found a 
number of new and unexpected features for the MMIT. 

(i) A local compressibility that appropriately char- 
acterizes Mott-insulating regions, shows critical behav- 
ior on entering those regions. Due to the microscopic 
nature of the phases, spatial correlations appear not to 
contribute to the critical behavior discussed here. This 
is a new form of MMIT, not observed so far in simple pe- 
riodic systems, that might be realized in fermionic gases 
trapped on optical lattices. Therefore, our observation 
adds a new aspect to this long-standing problem in con- 
densed matter physics. We expect that a similar local 
critical behavior will arise in higher dimensions as long 
as the spatial extent of the Mott domain remains finite. 

(ii) Universal behavior is found for k' for n ^ 1, in- 
dependent of the confining potential and/or strength of 
the interaction, excluding, however, the unconfined case. 
Also universal behavior is found for the variance A when 
n ^ 1. In this case, this behavior is shared by the un- 
confined model, (in) Finally, a proper scaling form for a 
characteristic density is introduced that leads to a generic 
phase diagram, with interesting features, as described in 
the paragraph above. 
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